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ABSTRACT 


A computer program was developed to obtain the wind 
tunnel wall corrections for wing angle of attack, induced 
Hl drag, and pitching moment in incompressible flow. The 
vortex lattice method is used for computation of these 
correction factors. The program can be applied to wind 
tunnels of arbitrary cross-sectional shape, and wings of 


any desired planform, subject to the constraint of straight 


leading and trailing edges. 


TABLE OF CONTENTS 


Ee INTRODUCTION -------~-~-------~-------~---~------ 
It. METHODS USED ---------~-------~-------~--~------ 
A. GENERAL --------~-~------------~---------- 
B. VORTEX LATTICE METHOD --------~-~----~----- 
C. AIRFLOW REPRESENTATION -------~-~---------- 


D. DETERMINATION OF WALL CORRECTIONS 
III. RESULTS AND DISCUSSION ----------~-~-----~------ 
A. THE COMPUTER PROGRAM --------~------------ 
B. GROUND EFFECT ~-------------~-~------~------ 
C. WIND TUNNEL APPLICATION --~-~--------------- 
Iv. CONCLUSIONS AND RECOMMENDATIONS ------------~- 
APPENDIX A: ASSUMPTIONS -----------~-~--------------- 


APPENDIX B: DEVELOPMENT OF INDUCED 
é VELOCITY EQUATIONS ---~--------------~- 


APPENDIX C: COMPUTER PROGRAM INPUT --------------~- 
COMPUTER PROGRAM -~--------------~--~----------~-~-~- 
LIST OF REFERENCES -<---------------------------~-~+--- 


INITIAL DISTRIBUTION LIST -------~--------------~--- 


et tee a ae ee ee, Se 


LIST OF FIGURES 


Wind Tunnel Geometry ---------~------~--------- 
Wing Geometry ------9-- 3-9 en ee = 


Ground Effect Ratio Plotted as a Function 
of Wing Height Above the Ground -----~--------- 


Wing Span Load Distribution --~------~--------- 
Induced Drag Distribution ----~------~--------- 


Angle of Attack Correction Factor Plotted 
as a Function of Span Ratio --~------~--------- 


Angle of Attack Correction Factor Plotted 
as a Function of Sweep Angle -~------~---------- 


Velocity Induced by a Vortex Segment ---------- 


18 


22 


32 
35 
39 


40 


41 
46 


LIST OF TABLES 


Ground Effect Factor Variation With 
Sweep Angle 


Wing Parameter Variation with Span 


Wing Parameter Variation with 
Sweep Angle 


[AERO] 


NS 


LIST OF SYMBOLS 
Wing aerodynamic influence coefficient matrix 
at the 0.75 chord position 
Wing aspect ratio = b/s 
Wing Span 
Wind tunnel cross-sectional area 
Local wing chord, streamwise direction 
Average wing chord, (s/b) 
Wing span loading, (%/q) 
Wing lift coefficient = Lift/qS, dimensionless 


Ground effect ratio 


itching moment coefficient = Moment/qSc 
»t chord 
Tip chord 


Wing induced drag coefficient 


Spanwise induced drag distribution (d/q) 
Wing induced drag distribution, drag/unit span 


Wing aerodynamic influence coefficient matrix 
at the wing quarter-chord 


Wall vortex loop strength 

Wing height above the ground 

Wind tunnel height 

Wing running span load, lift/unit span 


Number of horseshoe vortices per semispan 


Dynamic pressure 


Vet, 


Subscripts: 


FA 


Wall aerodynamic influence coefficient matrix 


Influence coefficient matrix for wall-on-tail 
Influence coefficient matrix for wall on wing 
Wing area 

Wing tape ratio, Cr/Ce 

Freestream velocity 

Induced downwash velocity 

Wind tunnel width 

Induced downwash angle at 0.75 chord position 
Influence coefficient matrix for wing-on-wall 
Influence coefficient matrix for wing-on-tail 
Axial coordinate 

Vertical coordinate 

Spanwise coordinate 

Angle of attack 

Wind tunnel wall interference factor 

Vortex strength 

Induced downwash angle at the tail 

Wing spanwise stations (33) 


Tail efficiency factor 


free air 
ground 


wall vortex loops 


wind tunnel 


a nner ong 


TL tail 
Ww wing 
( ie a9( )/da 


Matrix Notation: 


fi Rectangular matrix 
t 9 Column matrix 
[ee a Diagonal matrix 


[ ] Inverse of a square matrix 


ACKNOWLEDGMENT 


My sincere gratitude goes to Professor L. V. Schmidt 


: 
{ 


for his invaluable assistance throughout the course of this 


work and to my wife Ann2 for her patience and understanding. 


I. INTRODUCTION 


In the process of wind-tunnel testing the experimental 
procedure itself constrains and alters the flow past the 


wing under examination. This modified flow induces errors 


ithe inte i 


which must be corrected if accurate results are to be 
obtained. The major interference effects in subsonic wind 
tunnels result from the constraining influence of wind 
tunnel walls, disturbances from model supports and measuring 
equipment, and turbulence and non-uniformity in the airflow 
itself. 

Russell [Ref. 10] treats the problems resulting from 
model support interference and introduces a method for 
reduction of this effect. The differential equations 
governing flow within a wind tunnel are the same as those 
in free air. However, the presence of the wind tunnel walls 
changes the outer boundary conditions. The result is a 
series of perturbations of the velocity-potential gradient 
normal to the model surface. The effect is pervasive in 
that almost all of the measured aerodynamic quantities must 
be corrected. 

Generally, Prandtl is acknowledged as having established 
the foundations for research on wind-tunnel wall interfer- 
ence. His lifting-line theory and concept of trailing vor- 


tices opened the doors to much experimental investigation 


in this area. The resulting studies considered two-dimensional 


eset alts 


and three-dimensional effects, open and closed wind tunnels, 
and various geometrically shaped test sections. 

References 1 and 8 offer a comprehensive treatment of 
wind-tunnel wall interferences. Much tabular and graphical 
data are presented for different tunnel shapes and wing 
spans. Pankhurst and Holder included much data on octagonal 
wind tunnels, while Pope concentrated more on a uniformly 
loaded wing. 

With the advent of sweptback, slender wings designed 
for high-speed flight, the lifting line model was no longer 
suitable. In these situations, the wing lifting system 
may be represented by distributed horseshoe vortex elements 
(Refs. 6 and 7]. In Ref. 9, data, graphs, and formulas are 
easembled from various sources and systematically analyzed 
for various ude GE wings in wind tunnels. Although the 
mathematical formulations are tedious, the resulting graphs 
and tables summarized very well the data obtained from each 
theoretical treatment. 

Even with all this information available, it is still 
an arduous task to approach tabulated data with a specific 
wind tunnel shape and geometric configuration, and hope for 
a solution to the wall interference problem. It is the 
purpose of this paper to develop a computerized, closed- 
form solution to this problem. The design goal was to 
develop a digital computer code that would be applicable to 


wind tunnels of arbitrary cross-section and wings with 
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varying planforms. Such versatility would enhance the 


accuracy of test data and provide a valuable, accurate 


tool in wind-tunnel testing. 


II. METHODS USED 


A. GENERAL 
Wind tunnels provide an invaluable tool in the design 
and performance prediction of wings and wings-and-bodies. 
However, the wind-tunnel testing environment introduces 
its own set of problems. The interference action from the 
walls of the wind tunnel alters the development of induced 
velocities by the winglifting system. As a result, downwash 
formation is suppressed and the wing aerodynamic character- 
istics are those of a wing having greater aspect ratio. The 
accurate determination of wind tunnel wall interference is, 
therefore, an important consideration in wing experiments. 
The classical wall correction factor, 5, for the case 
of a single horseshoe vortex representing the complete wing 


is defined by: 


is 
Q 
" 
OH 

aQIn 
Q 


(1) 


where 


Using small angle approximations, the tangent is approxi- 
mately equal to the angle. The wall-correction factor may 
be expressed in terms of wing circulation, vortex span and 


average induced downwash velocity by: 


15 


S$ = 25 (2) 

Calculation of 6 is dependent on the determination of 
the downwash velocity induced by the wind tunnel walls. 
Classical methods involve the utilization of an image 
system of vortices which cancel the effect of the wing 
lifting system at the wall. The boundary condition of no 
flow through the wall is, therefore, preserved. This system 
can provide an accurate, closed-form solution to the ground 
effect problem. However, when four walis are present, each 
wall must be represented by an image system which affects 
other image systems, as well as the model itself. The 
result is a series of infinite systems. Truncation of this 
system introduces inaccuracies of unknown magnitude and 
significance. 

This method of images has been applied to tunnels of 
rectangular cross-section with success. However, it is 
inadequate when arbitrary tunnel cross-sections are 
encountered. Add to this a wing with sweep and taper, and 
the problem becomes significantly altered. 

Therefore, a new method was sought to solve this problem 
specifically for the NPS 3.5 X 5 ft. octagonal wind tunnel. 
The technique chosen was presented by Joppa in Refs. 4 and 
5. This approach lends itself very well to computer solution 
and can be applied to wind tunnels of any cross-sectional 


shape. 
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In this study, an extension of Joppa's vortex lattice 
method was made to allow the treatment of wings with varying 
shape, sweep, and aspect ratio. Such treatment requires a 
modification of the single horseshoe vortex representation 


of the wing lifting system. The method of Weissinger as 


aL ee 


presented in Ref. 6 provides an accurate means of representing 


variable planforms. 


ee ie 


B. VORTEX LATTICE METHOD 

Since the method of images was considered to be limited 
in application, it was necessary to find another, more versa- 
tile method of representing the wind-tunnel walls. The 
problem is essentially one of finding a vorticity distribu- 
tion which will accurately model the wind-tunnel walls. This 
vortex system must satisfy the boundary condition of no 

] flow through the walls. 

In Ref. 4, Joppa presents a method which accurately 
models the wind-tunnel walls. This was accomplished by 
representing the tunnel walls with a lattice network of 
rectangular vortex loops. This vortex system lies in the 
plane of the tunnel walls and satisfies Helmholtz's theorem 
that a vortex filament can neither begin nor end in the 
flow. Fig. 1 shows a representation of this system as 
applied to the NPS 3.5 by 5 ft. wind tunnel. 

Each vortex segment in a given rectangle has vortex 


strength Tie and has a corresponding control point located 
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in the center of each rectangle where the no-flow boundary 


condition is satisfied. The last vortex lattice rectangles 
located far downstream are represented by horseshoe vortices 
trailing to infinity. The reasoning is that far downstream, 
where the influence upon the wall vortex lattice array of 

the bound vortex elements on the wing is small relative to 
the trailing vortex elements, one finds that the vortex 
strengths of the lattice terms do not change in the stream- 
wise direction. Consequently, by a straightforward cancella~ 
tion argument, one can visualize replacing the vortex lattice 
rectangles by horseshoe vortices trailing to infinity. 

The mathematical presentation of the theory is based on 
the Biot~Savart law describing the velocity induced at a 
point by a vortex segment. Joppa in Refs. 4 and 5 follows 
the mathematical formulation beginning with the Biot~Savart 
law. This development is summarized in Appendix B. 

The general matrix representation of this method centers 
around the boundary condition of no flow through the wall. 
As mentioned previously, this condition is satisfied at a 
control point located at the center of each vortex loop. 
Therefore, it follows that the velocity induced at each 
control point by the vortex loops must be exactly cancelled 
by the velocity created by the wing lifting system. This 


is presented below in matrix representation 


(ROL] {Pr} = (WOL] ity! (3) 


which allows a solution of vortex lattice strengths in 


terms of the wing vortex strengths. 


_ =1 
{Th} = [ROL] 


[WOL] tht (4) 
Note that the [ROL] matrix, although possibly quite large, 
is square and hence may be inverted. 

{ROL] represents the wall influence coefficient matrix. 
It is obtained by calculating the downwash induced at each 
wall control point by a unit value of each vortex loop. 
Therefore, Bele is the downwash at wall control point "i" 
induced by vortex loop "j". [WOL] is the influence coeffi- 
cient matrix describing the wing effect on each loop con- 
trol point. Wok is the downwash produced at control point 


"i" by the — 


wing horseshoe vortex. 

By inverting [ROL] and premultiplying with [WOL], the 
strength of the wall vortex loops can be found in terms 
of the wing circulation. The inversion of [ROL] can be 
very time-consuming, since it may bea very large matrix. 
However, by using symmetry, the order of this matrix can 
be reduced by a factor of four. This requires the wing to 


be placed on the horizontal and vertical centerline of the 


wind tunnel. 


C. AIRFOIL REPRESENTATION 


In Ref. 5, Joppa stated that the wing can be repre- 


sented by a lifting system more complicated than a single 


elementary horseshoe vortex. In order to model swept wings 
of arbitrary planform, a more definitive lifting system 
model was employed. Gray and Schenk, Ref. 6, present a 
modified Weissinger method designed to find the load dis- 
tribution of wings with arbitrary planform. The wing is 
represented by a series of horseshoe vortices. The midspan 
of each horseshoe vortex bound element is located at the 
wing quarter-chord position. The condition of flow tangency 
is satisfied at the wing three-quarter-chord point in the 
center of each horseshoe vortex. Figure 2 illustrates the 
configuration chosen. 

Pearson [Ref. 7] modified this technique by allowing 
the bound vortices to be swept, and including chord-wise 
horseshoe vortices. For this presentation it was felt that 
unswept bound vortices would accurately represent the lifting 
system of a swept wing. In addition it was felt that little 
deterioration of accuracy would result if the chordwise 
system of bound vortices was combined into a single bound 
vortex at the quarter-chord. Figure 2 illustrates the 
configuration chosen. | 

An important concept in this representation of the wing 
is that of flow tangency at the three quarter-chord. The 
downwash angle at the wing can, therefore, be determined 


by placing a control point at this position in the center 


of each horseshoe vortex. The rationale behind the three- 


quarter chord flow tangency condition is a product of 


pth ce 


Figure 2. Wing Geometry 
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two-dimensional flow considerations. Extension of this 
technique to the three-dimensional wing has been shown to 
produce experimentally confirmed results. However, it is 
still unclear why this technique works so well. 

The wing aerodynamic influence coefficient matrix was 
then developed. [AERO] represents the downwash induced at 
each 0.75 chord wing control point by a unit value of wing 
horseshoe vortices. The matrix formulation for induced 


downwash formation in free air is presented below: 


{w} = [AERO] Ces, (5) 


Dividing the above expression by freestream velocity, one 
obtains an expression for 0.75 chord control point induced 
downwash angles. The flow tangency requirement equates 


this column matrix to the wing angle of attack column 


vector. 
8} = [AERO] {5} (6) 
3c/4 
and 
fa} = {3} 
3c/4 
It should be noted here that < is equal to from 


the Kutta~Joukowski Law. Substitution into equation 6 


yields the following results: 


z3 


(@) = 0.5 tazRo) (4) (7) 
i QW pa | 


Solving equation 7 for (Sh, one can obtain the wing span 
load distribution. Since = is equivalent to the sectional 
value of Coc, integration oe the = distribution vs. span 
provides a direct measure of the wing lift coefficient, Che 
When one assumes a uniform angle of attack value of one 
radian along the span, the resultant integration of the 
span load solution will yield the wing lift curve slope, 
Cra? It should be borne in mind that these techniques are 
predicated upon principles from thin airfoil theory and 
in the limit for an infinite aspect ratio wing will yield 


ul 


a wing lift curve slope of 2m(rad ~). 


D. DETERMINATION OF THE WALL CORRECTIONS 

The classical wall interference factor 6, which includes 
a ratio of wing to tunnel cross-section area, is generally 
found by representing the wing as a single horseshoe vortex. 
Since the present formulation utilizes ten horseshoe vor- 
tices per semi-span, § in a modification can best be found 


by continuing our matrix algebra. 


By assuming that the lift coefficient is the same in 


both tunnel and free air, the following equations result: 


+ 6 cL (9) 


(10) 


The free air lift slope was found in the previous section. 
The wing lift curve slope for the wing in the wind tunnel 
can be found by a simple extension of this same method. 

The total downwash at the wing control points is equal to 
that induced by the wing added to that created by the tunnel 
walls. The effect produced by the tunnel walls is repre- 
sented by the vortex rectangle-on-wing influence coefficient 


matrix [ROW]. iar is the downwash induced at the ee wing 


control point by the — 


wall vortex loop. Multiplying this 
matrix by the vortex loop strength, one obtains the total 


wall effect. 


{wi} = [AERO] uF + [ROW] (f,} 


Substitution of equation 4 for {r,} forms the following 


set of equations: 


a 


{wi} = [AERO] (To + [ROW] [ROL] ~ [WOL] iT 3 


1 


(3) = 0.5 [[AERO] + [ROW] [ROL] two] ] {4} (12) 
I 


By establishing one radian angle of attack, solving 
for = and integrating the result, the wind-tunnel lift 
curve slope can be determined. It should be noted that the 
correction factor obtained from this process is different 
from the classical $5, since the correction factor computed 
above already contains the wing area to tunnel area ratio 
influence. Therefore, the correction to angle of attack 


can be obtained by a direct application of this 6. 
Nas =_ 656 (13) 


In most literature, the correction for induced drag 


is obtained by application of the same 6 as that used for 


angle of attack. 


However, in this presentation, the induced drag correc- 
tion is calculated in a different manner because of added 
available aerodynamic information. The spanwise induced 
drag distribution in free air is given by the following 


formula: 


P 
d = (<4) [DAERO] {2} (15) 
q Pa qq” FA q aM 


{[DAERO] is the wing aerodynamic influence coefficient 
Matrix computed with the wing control point now located 
at the quarter-chord position. The wall effect can similarly 


be found by including additional matrix terms for the wall 


effect on the wing quarter-chord position. 


(3; ee [[DAERO] + [ROW] [ROL] ~* {won} 1 {2} (16) 
q q Ip 


BY on 


Integration of equations 15 and 16 result in the free 
air and wind tunnel induced drag, respectively. The induced 
drag correction factor can then be found in the following 


manner. 


2 
“Dy ~ Bern “legs 
FA 
7 2 
“D, oy a 
T 


are aise ei AEF 


colt EE 2 
Kea = py. foes 
tPA 
ce 2 
Kp = Cy /Cy, 
lop at 
(iz) 
K = Rea Re 
2 
AC. = KC 
D; Lin 


When testing a model in the tail-on configuration, 
the presence of the wind tunnel walls serves to modify the 
downwash angle induced at the tail. This altered flow 
field can have a significant effect on the pitching moment 
coefficient. The following increment in pitching moment 


coefficient is generated by the horizontal tail: 


AC se, eee LS (18) 


The term in parentheses is approximately equal to the 
3c 

stabilizer effectiveness, 35" . The tail angle of attack 
i 


can be estimated as follows: 
a = @= ¢ + 6. (19) 
si 


At a given value of wing lift coefficient an incremental 
change in downwash angle at the tail due to the wall effect 


can be computed: 
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Ae = €. =e (20) 


Therefore, the correction to pitching moment coefficient is 


given by the following relation: 
Ac = €.. (21) 


Since the downwash angle €« is equal to 7 at the tail 
position, e« for free air and for the wind tunnel can be 


found from the following relations, respectively: 


Q 
= 0.5[(woT) {= 
Cpa [ ] eo 


e, = 0.5{{WOT] + [ROT] [GAMAL]] oe (22) 
[WOT] is the wing effect on the tail position and 
{ROT] [GAMAL] represents the wall effect on the tail. 
Equation 20 can then be used to find Ae. The pitching 
moment correction, as evaluated in equation 21, is thus 
dependent on being given an experimental value for the 


tail effectiveness parameter. 
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III. RESULTS AND DISCUSSION 


A. COMPUTER PROGRAM 

The computer program developed consists of two major 
subroutines. Subroutine Loop calculates the downwash at 
any point due to a vortex rectangle. Subroutine HSHOE 
computes the downwash at any arbitrary point due to a 
horseshoe vortex. HSHOE is used for both the wing horseshoe 
and those of the wall system located far downstream in the 
tunnel walls. The program offers two basic options. One 
option establishes whether ground effect or wind-tunnel 
wall corrections are to be calculated. The other option 
indicates whether the tail-on or tail-off case is to be 
considered. 

Each computer run offers the inclusion of many wings 
of different span, aspect ratio, taper ratio, and sweep. 
The wing parameters are read as input data after the inver-~ 
sion of the large wall-on-wall coefficient matrix. There- 
fore, with a given tunnel configuration, a variety of 
different wings can be included with a minimum of computer 
time consumed. | 

Integration of the wing span loads in accomplished by 
a subroutine obtained from the NPS source library. It 
should be noted that in the integration scheme the assumption 


is made that the wing loading at n = 0.0 is the same as that 


for n = .05b. Any inaccuracy resulting from this assumption 
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appears to be negligible. Appendix C describes the necessary 


input to the program. 


B. GROUND EFFECT 

Originally, ground effect calculations were conducted 
merely to check the accuracy of the vortex lattice method. 
As previously stated, the method of images provides a 
closed-form solution to the ground effect problem. There- 
fore, this provided an excellent opportunity for a conver- 
gence check. As an object of academic interest an option 
for ground effect computation was retained in the final 
conputer program. 

The results of the convergence check showed that the 
vortex lattice method could very accurately compute the 
velocity induced at a point by the presence of a wall. By 
applying the theory previously developed, the lift curve 
slopes for free air and ground effect were computed. Results 
were then compared with those presented in Etkin [Ref. 3]. 
The free air lift curve slope was found to correlate extremely 
well for both straight and sweptback wings. 

Furthermore, in Ref. 3, Etkin defines a ground effect 
ratio as the ground effect lift curve slope divided by that 
in free air. The results are tabulated in Table 1, and 
Figure 3 compares the results obtained from the computer 
program with those presented in Etkin. As can be seen, 


agreement is very favorable. 
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Figure 3. Ground Effect Ratio Plotted as a Function 
of Wing Height Above the Ground 
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Several computer runs were conducted at different wing 
heights above the ground. For each run, sweep was varied 
from 0 to 40 degrees. Results are tabulated in Figure 4. 
It can be noted that the ground effect ratio initially 
increases with sweepback but returns to the unswept value 
at about 30 degrees. 

As a note of caution, care should be exercised when 
the wing gets very close to the ground. It is necessary 
to greatly decrease the vortex loop size for accurate 
results. This procedure seems to provide a much more finely 
defined ground representation. Otherwise, calculations 


will result in inordinately high values. 


C. WIND TUNNEL APPLICATION 

A convergence check of the wind tunnel wall correction 
factor was obtained by applying the program to a 3 X 5 ft 
rectangular tunnel. The correction factor obtained was 
converted to the classical 6 by including the factors 2 
and 57.3 (radians to degrees). The result of 0.1216 was 
compared to that presented by Joppa in Ref. 5. The wing 
geometric configuration presented in Ref. 5 consisted of a 
single horseshoe vortex. The value of 0.127 from Ref. 5 


correlates very well with the results of this study. 


Tables 2 and 3 present the data obtained from compu- 


ii lta a 


tations with the NPS 3.5 X 5.0 ft octagonal wind tunnel. 


Computer runs were conducted using wings of varying span 
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Figure 4. Wing Span Load Distribution 
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and sweep angle. Figures 4 and 5 illustrate the span loads 


calculated for a typical wing. It is interesting to note 
that as sweep angle increased, negative induced drag appears 
on the outboard wing stations. Therefore, total induced 
drag is seen to steadily decrease with increasing sweep 
angle. 

As sweep angle is increased, the free air lift curve 
slope decreases. Correspondingly, the corrections to 
angle of attack and induced drag are reduced. Figure 6 
graphically displays the data for a straight wing. In the 
past, the correction factors for the NPS 3.5 X 5.0 ft 
octagonal wind tunnel were obtained by assuming an elliptic 
cross~section. The correction factors presented in Ref. l 
for elliptic cross-section are 0.610 for angle of attack 
and 0.0106 for induced drag. In comparison, the vortex 
lattice program computes values of 0.7086 and 0.01075, 
respectively. It is interesting to note that the induced 
drag correction is almost identical, whereas the correction 
to angle of attack is slightly higher for the octagonal 
wind tunnel. Figure 7 illustrates the angle of attack 
correction factor with varying tunnel span ratios. The 
correction factor can be seen to steadily increase with the 
span ratio. 

It is an interesting observation that the value of 


induced drag for a wing span of 5.0 ft was almost zero. 


This compares very well with airfoil theory in that a wing 
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Figure 5. 
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Figure 6. Angle of Attack Correction Factor Plotted 
as a Function of Span Ratio 
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Figure 7. Angle of Attack Correction Factor Plotted 
As A Function Of Sweep Angle 
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which spans the wind tunnel represents an approximation to 

the two-dimensional wing. Induced drag for a two-dimensional 
wing should equal zero. In addition, the wind-tunnel lift 
curve slope is within 0.0006 of the theoretical two-dimensional 
value of 27 raat, 

A sample computer run was made for the moment correction 
factor. The tail quarter-chord was placed two feet behind 
the wing bound vortex, on the tunnel centerline. The 
correction factor obtained for a straight wing of aspect 
ratio 6.0 was -0.2297. For a wing with sweep angle of 
30 degrees and taper of 0.3, a value of -0.1887 was calcu- 
lated. The above correction factors must be multiplied 
by the stabilizer effectiveness parameter, using consistent 
angular measures, in order to find the pitching moment 


correction. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 


It is felt that the vortex lattice computer program 
as presented in this paper offers a very reliable and 
flexible tool for the determination of wind-tunnel wall 
corrections. It is capable of treating wings of any span, 
taper ratio, aspect ratio, and quarter-chord sweep, subject 
to the constraint of straight leading and trailing edges. 
The geometric shape of the wind tunnel cross-section can 
Similarly be varied at the discretion of the program user. 
Resultant accuracy appears to be excellent. However, cau- 
tion should be exercised in determining the vortex loop size. 


As the wing is positioned closer to the walls, a smaller, 


‘ more definitive vortex loop representation must be used. 


A logical extension of this program would be the inclu- 
sion of wings with an antisymmetric load distribution. 
Furthermore, slight modifications could incorporate the 
ability to consider the wall effects on other stability 
derivatives. In its present form, this program offers 
the wind tunnel user a versatile tool for determination of 


wind tunnel wall corrections. 
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APPENDIX A 


The general assumptions used in the compilation of 


this program are listed below. 


E. 


2. 


The two-dimensional lift-curve slope is 27 raat, 


The wing load distribution is symmetrical. 
Aeroelastic effects were not considered. 

The wing is located on the horizontal and vertical 
centerline of the wind tunnel. 

The leading and trailing edges of the wing are 


straight. 
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APPENDIX B 


DEVELOPMENT OF INDUCED VELOCITY EQUATIONS 


The equations required to solve for induced velocities 
consist of two main types. One equation solves for the 
velocity induced by an arbitrarily oriented vortex segment, 
whereas the other basic equation applies to horseshoe vortex 
elements. These equations have been incorporated in sub- 
routines LOOP and HSHOE, respectively. The development of 
these equations presented below follows closely that found 
in Ref. 4 and 5. 
According to the Biot-Savart law, the velocity induced 


at a point p is found from: 


eae worn 


Us 
aah (cos B 


=| 
t 


l + cos Bo) Vv (1B) 


For vortex segments oriented in the Y-Z plane, 


Be ok 


Mie: Stet: ee er ee 
cos By + cos B, = 7R, R55 (S (Ry R,) ) (2B) 


It is now necessary to consider v, the unit vector 
which establishes direction. From Fig. 8, it can be seen 


that 


Ft (3B) 


! | 
; 
. E| 
; bi 
' 
(XPYPZP) 
| Figure 8. Velocity Induced By A Vortex Segment | 
] 
i 


Noting from Fig. 8 that 


hs = |R, x S| (4B) 
then 
he eo ee SS Se ee, ce 
ly z lz y 1z°x LXSZ 
: 2, 1/2 
+ (Ry 8, Ry ySx) ] 


Calculating the induced velocity in the direction of 


outward normal of point P per unit I, 
; = = -n (5B) 
where the unit normal vector is defined by: 
n = oi + (sin 9p)j + (cos 9p)k 
Therefore, one finally obtains: 
R, + Ro 2 


= ——+—_4, Is* - (Ry ~ R3)7] 
87R, R, (hs) 


j= 


(6B) 


* [(sin Op) (R, 2S, - R),.5,) + (cos ep) (RS - R/S.) ] 
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A conditional statement must be included in the program 


to set the downwash equal to zero if hs is equal to Zero. 


This will preclude having a zero term in the denominator. 


The trailing elements of a vortex loop can be treated 


with more simplicity. For these segments, 


(R,? = n2)2/2 
cos B 
1 Ry 
(7B) i 
(R,? 2)1/2 
; cos B = 
j 2 Ry 
i However, it can be seen from Fig. 8 that 
i" = (xp - xA)? + (yA - yp)? + (2a - zp)? 
(8B) 
hn? = (yA - yp)? + (2a = zp)? 
Combining equations 7 and 8, 
cos By = xp 
1 
(9B) 
cos Bo = AOE 
2 


As a result, the equation for the trailing elements 


of a vortex loop is simply 


cos B, + cos B 


= —1,—_ {(sin 6p) (zA - zp) 
4th 


(108) 


+ (cos 9p) (yp - yA) ] 


For the horseshoe vortex, the equation for the bound 
segment is identical to equation 6. Similarly, the 
trailing elements can be treated as in equation 1l with a | 
slight modification. In this case, the cosine of one 
encloseu angle is equal to 1, since the trailing element 
extends downstream to infinity. The induced velocity 


equation becomes 


| w 1 + cos By wa - 
| Foy 2A - 2p)5 + (yp - yA)k) 
4th 
| 
Ww 1 + cos B 

T = 


——,—I (sin @p) (ZA - zp) 
4th 
(11B) 


+ (cos Op) (yp - yA) ] 


Card Number 
1 


Input cards 5 and 6 can be repeated for different 


Pplanforms and tail locations. 
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APPENDIX C 


COMPUTER PROGRAM INPUT 


Column Number 


1-72 


Description 


Program title 


MM=Number of Spanwise 


Stations 


NN=Number of Downstream 


Stations 
Deltax=Loop length 


NOPT=l1 for Ground 
Effect 


NTAIL=1 for tail-on 


Tunnel y-coordinates 
(one guadrant only) 


Tunnel z-coordinates 
(one quadrant only) 


Span 

Aspect Ratio 

Taper Ratio 

c/4 Sweep (degrees) 


x-Coordinate of tail 
quarter-chord 


z-Coordinate of tail 
quarter-chord 


Format 
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